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SUMMARY 


A study was undertaken to examine present and proposed methods of 
calculating the attenuation afforded "in and down" scattered radiation by a horizontal 
slab, and the effects of finite fields of contamination in order to recommend improve¬ 
ments and up dating for OCD documents TR-20, Volumes I and 2 , "Shelter Design 
and Analysis." 

In this study we examined the assumption used in computing the "in and 
down" attenuation factor, assessed their effect on the resulting attenuation factors, 
and compared the results with the latest available experimental data. Recommendations 
are provided as to "best value" factors and how they might be included in existing 
publications. 

The effects of finite fields of contamination were subjected to further 
analysis by directly comparing the existing method of calculation with results 
obtained using transmission coefficients calculated by the Monte Carlo method. This 
comparison, though complete only for above-ground locations, indicates that further 
study is required before the existing method of calculation can be relied upon. 

The study of in and down scattered radiation is described in Chapters 1 
through 3: finite fields of contamination are in Chapter 4. These major conclusions 
may be drcwn from this study: 

The "In and Down" Attenuation Factor 


1. The in and down attenuation factor, (B' 0 (X, u), represents a significantly 
better estimate of this effect than the values previously used (Bg (X) ). 

2. Calculated and measured attenuation factors agree with : n 18.5 per cent; 
the analytical value is conservative. It was found that the standard 
deviation of the ratio of the attenuation calculated, Bg(X,u), to that 
measured in 346 experimental measurements, was 13.5 per cent. 


Finite Fields of Contamination 


3. The current method of computing the effects of finite fields of 
contamination requires further study. 

20 

In total dose, the NBS Monte Carlo compilation method and 
that presently proposed ^ agree fairly well for detector locations 
on the first floor of a sheltering structure. 


4. 
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5. Major differences in the fraction of scatter to direct dose were 
found when using the NBS method, casting serious doubt on the 
accuracy of current methods of shelter analysis when the detector 
is shielded from direct radiation (i.e., in a basementor an 
upper-story location shadowed by a heavy floor). 
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CHAPTER I 

IN AND DOWN SCATTERED RADIATION - THEORY 


INTRODUCTION 

Most of the error observed between measured and computed values of 
dose :ote »n basement areas has been attributed to the attenuation provided by the 
basement ceiling. Values of this attenuation previously were computed in approxi¬ 
mate form for cobalt, cerium, and the 1.12 hr. fallout energy spectrum, using a 
single Compton scatter energy change and transmission coefficients computed by 
the Monte Carlo method, assuming that all incident radiation lay in the plane 
mutually perpendicular to the exterior wall and the ground and that it passed through 
the detector location. Agreement between the attenuation computed in this 
approximate fashion and that measured was much better than that obtained with the 
engineering method. 

Our purpose in this study has been to re-examine the effects of the 
assumptions used and to recalculate the basement ceiling attenuation more accurately 
by (1) relaxing the restriction that incident radiation must lie in the plane containing 
the detector that is mutually perpendicular to the ground and the wall and (2) using 
a finer gradient when performing the numerical integration over the contaminated area. 


Basement Ceiling Attenuation 

Approximately 80 per cent of the radiation arriving at the exterior wall 
of a structure located in an infinite contaminated field is uncollided. The remaining 
20 per cent is singly or multiple scattered radiation of lower energy. This "atmospheric" 
scattered radiation is divided about evenly between the upper and lower hemisphere 
for a detector located above ground level. Thus, in estimating the attenuation afforded 
by a basement ceilings the first assumption is that all radiation arriving at the 
above-ground v-alis of a structure that scatters to the basement area is uncollided. 
According to this assumption, all scattered radiation is of higher energy, therefore 
rttenuating less rapidly tnan actual scattered radiation, hence producing a conservative 
estimate of the radiation scattered to a below-ground area. To simplify the problem, 
we adopted a second assumption: radiation scattered by the above-ground portion 
of the structure is reduced in energy by a single Compton scatter, through the angle 
required to reach the detector only. This assumption is also conservative in that it 
treats all multiply scattered photons as of higher than their actual energy. Raso has 
computed the average energy of scattered photons for incident energy of 1.25 mev 
and finds that this figure is about 85 per cent of that predicted by single-scotter 
theory. The third major assumption is that the detector is relatively near the basement 
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ceiling; i.e., radiation arriving at a basement detector by a ceiling slab is attenuated 
exactly as much as that of parallel monoenergetic radiation of the same energy and 
angular incidence in relation to the external surface of the basement ceiling. 

The procedure used to compute the ceiling attenuation was as follows 
(see Figure 1). The contaminated field was subdivided into "n" increments of equal 
radiation worth, for each of the nine energy components of the 1.12 hr. fallout 
spectrum, when viewed by a detector located at the scattering area. The total angle 
of scatter (\jj) then was determined at the scattering location between the source 
patch and the detector. Notice that the previously used model for computation 
was a cylinder (or the equivalent, the short wall of a building with large eccentricity) 
and therefore the exit azimuthal angle (3 was identically zero. Then, the new 
scattered energy was computed from the Compton scatter equations. The dose 
arriving at the detector location is taken as proportional to the amount of radiation 
energy arriving at the wall from the source patch, reduced by the ratio of the initial 
energy to that found after scattering, and multiplied by the attenuation introduced 
by the basement ceiling (taken as equal to that of parallel monoenergetic radiation 
of identical energy and angular incidence) and the air absorption coefficient. 

A quantity proportional to the total overall dose may then be calculated by summing 
this quantity over all energies of the fallout spectrum weighed by the energy spectrum 
and over the infinite plane contaminated area. The ceiling attenuation factor was 
then computed by dividing this quantity by the similar quantity for zero ceiling 
thickness. 


This may be expressed in equations in the form presented here (see 
Figure 1). The contaminated field beyond the wall is subdivided into a total of "n" 
contaminated patches of equal worth (when viewed by a detector located at the wall 
scattering volume) in cylindrical geometry. The total uncollided dose rate at height h 
from a semi-infinite field of contamination is: 


D(E) = const 


.- |i t (E)p 2,pd p 


const Ej(p(E)h) 


.where: 


M t (E) 


P 

E, 


= is the total coefficient in air for photon energy E 
= the height of the scattering element 
= the slant distance from detector to differential source area 
= the exponent integral of the first kind 
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The total number of differential patches "n" is composed of "i" rings of 
equal worth, each subdivided into "j" azimuthal patches. The center of each patch 
thus is specified by coordinates 9| and a^ relative to the scattering area. 


a 

m 


= the dzimuthal coordinate (see Figure 1) of source patch "m" 


a 

m 


it(m + 1/2) 
1 



the polar coordinate of a source patch for contamination of 
ehergy E 


H t (E)h 

0j(E) = cos" 1 H t (Srp^ 


E^ t (E)P;) = 


E^M f (E)h) 


T 

I 


tt+ 1/2) 


For heights of the scattering volume "h" that are small with respect to 
a mean free path, the quantity 9.(E) is nearly independent of "h". 

I ' 

The detector is located with respect to the scattering volume by azimuthal 
angle (3 and polar angle 9. The cosine of the total angle through which a photon must 
scatter from the source patch "I", "m", to reach the detector therefore is: 


cos 4-' 


scos 9 cos 9. - sin 9 sin 9 ft cos (180 - a 

A { X m 



which is implicitly a function of the initial energy E 0 . Thus the energy E in mev of 
the radiation scattered b> the wall toward the detector, after suffering a single Compton 
scatter, would be 


E = E 

o 


.TRT 


.51 

TT 


cos 40 


;' 

•j 
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If the probability that a photon will scatter from source patch 9|,a^ to 
a detector located at 9, (3 relative to the scattering volume, is written as 
P(0|, Qm -'*■9,(3), the dose seen by such a detector located below a floor of thickness X 
from a single Compton scatter is 


D(E q ,X,9,(3) 


A£ p{9„ a 9, (3) E p n (E) cos 9 S'(X, cos 9,E) 
o l m u 

n E 

o 


where: 


AE = the total quantity of energy E q incident on the scattering 
volume 

E q = the energy of interest 

E = the energy of the photon after suffering a single scatter to the 

detector 

p^(E) = tissue equivalent energy absorption cross section 

cos 9 S'(X, cos 9, uE) = attenuation factor for parallel monoenergetic 

radiation of energy E, angle of incidence 9, 
and slab thickness X 

The attenuation introduced by the ceiling above the detector therefore 
is equal to the integral of this quantity over the total scattering volume (the wall) 
divided by the same quantity, with thickness X equal to zero. In practice the 
quantity P(9|, _»9,(3) is nearly constant and may be removed from under the 

integral sign and canceled. If P(C|, c^-^9,(3) is taken as approximately equal '■o 
the differential Compton scatter term rather than constant, the differences are quite 
small. This is Plustrated in Table 1 for the region of most interest, cos 9. Here 
agreement is found to be well within the region of other uncertainties of the calculation. 
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TABLE 1 


Comparison of attenuation factors computed with probability of scatter 

equal to the Compton differential scatter probability 

and a constant for (3 * 0 and E * 1.25 mev 

o 


cos 0 

1.0 

.75 

.50 

.25 

___ 

.10 

Scatter 

probability 

Comp. 

Const. 

jjjjjjjj 

Const. 






mu 

X f 

10 psf 

.90 

.89 

.83 

.82 

.68 

.67. 

.39 

.36 

.13 

.12 

20 psf 

.76 

.72 

.64 

.62 

.45 

.42 

.18 

.16 

.035 

.032 

40 psf 

.49 

.43 

.37 

.38 

.19 

.18 

.052 

.050 

.010 

.0091 

80 psf 

.18 

.12 

.11 

.029 

.039 

.038 

.0086 

.0070 

.0018 

Bb 

160 psf 

.023 

.010 

.0098 

.0068 

.018 

.0014 

.0038 

.00030 

.00081 



ror simplicity, the scattering probability may therefore be taken as constant 
over the range of interest of this problem, and the attenuation at position 0, may 
be written as 


Att (9, h, X, E o ) 


\ p (E) cos OS'fo cos 0,E) 
c a 
o 


| H (E) cos 9S'(X = 0, cos 0,E) 
c a 
o 


Two aspects of the problem may now be evaluated quite easily.,. First, 
since Raso^ found that the average energy of scattered photons was about 85 per cent 
of Compton single-scatter energy, what effect would this degradation of energy of 
the scattered photons have? Second, what effect other than the straightforward 
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geometric difference of slant penetration distance in the floor would a variation 
in angle beta have on the attenuation factor? 

These effects are illustrated in Table 2. Here again it is evident that 
essentially there is no secondary geometric effect on the floor attenuation and that 
the assumption of scattered energy equal to 85 per cent of Compton single scatter 
lowers the attenuation factors by only a small amount. 

Since most buildings are other than cylindrical in shape the necessary 
integration over the variable beta must be examined. Previously, it was assumed 
that the angular distribution of radiation emerging from the rear of a vertical wall 
in the plane parallel te> the ground is closely approximated by the cosine of the 
emergent angle relative to a perpendicular to the scattering surface. It was demonstrated 
above (see Table 2) that the principal effect of this angle being other than zero is 
in the straightforward geometric increase in the slant penetration distance of the 
slab over the detector that a photon must traverse to reach a detector a given distance 
from the wall. This distance, represented equivalently by an increase in the angle 9 
to 9 1 for an attenuating slab of constant thickness, may be calculated quickly by the 
relationshi p: 


Tan 9‘ = tan 9/cos |3 

The attenuation afforded by a rectangular structure may, then, be calculated in a 
manner similar to that applied to a cylindrical structure by integrating the dose 
scattered from a given scattering volume (represented by angle |3), weighed by the 
emergent angular distribution for a finite thickness slab divided by rhe similar quantity 
for a slab of zero thickness. Thus 


- 


all walls 

i 

i 


AD (9 1 , |3,X, E ) cos (3d (3 

o 


all walls 


AD (9’, (3, X = 0, E ) cos (3d (3 
o 


where 8^(X, Q, Ep, ^ ) - the attenuation that would be afforded by a slab to in and 
down scattered radiation in a rectangular structure. 
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L/W = length-to-width ratio of the structure 

This quantity has been calculated for several conditions of floor mass, 
building size, and eccentricity by numerically summing over ^he azimuthal dis¬ 
tribution in 5° increments. The results of this series of computation ure shown 
in Table 3, where they are compared with the values obtained by computing the 
floor attenuation of a cylindrical structure of equal solid angle fraction when 

TABLE 3 

Comparison of attenuation factors computed in circular 
and rectangular geometry for equal sol id angle fraction, u 
E = 1.25 mev 

(Note Structure is of differential wall height) 



viewed by the detector. This table indicates that these two methods give excellent 
agreement and that there is therefore no significant "geometry factor" for floor 
attenuation. 

The assumptions made in calculating the attenuation afforded to in and 
down scattered radiation by a floor slab upon evaluation, have been found to 
effect but Sidle the predicted attenuation values. Using a constant-scatter probabili 
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rather than the Compton single-scatter probability slightly decreases the predicted 
basement detector dose rates, whereas assuming scatter energy equal to the single 
Compton scatter energy and neglecting secondary geometric effects in the rectangular 
wall case slightly increases predicted dose rate. 

The floor barrier computation as described above has been performed 
for the approximate energy of cesium and cobalt radiation as well as for each of 
the nine energy components of the 1.12 hr. fallout spectrum. Figures 2 through 4 
present the results of this series of computations as a function of the mass thickness 
of the attenuating slab and the solid angle to the scattering volume, in Figure 4 
is presented the result of summing the data for the nine energy groups of the 1.12 hr. 
fallout spectrum weighed by the value of each of these groups. Previously estimated 
values produced by a rough hand calculation are generally within about 10 per cent 
of the values shown. 

The use of these figures is quite straightforward. In calculating the exact 
dose in the basement area, the dose from each differential increment of height of 
the scattering wall should be multiplied by the attenuation factor (Figures 2 through 
4) for the basement ceiling at the solid angle fraction of that differential wall height 
(when viewed by the detector) and these results summed over the entire scattering 
wall height. It was demonstrated earlier, however, that an equivalent value, 
generally within 10 per cent of that calculated by this method, may be obtained 
easily by averaging the solid angle to the top and bottom of the scattering wall and 
determining the attenuation factor for this average solid angle fraction. Several 
examples of this procedure (for a structure of 10 ft. wall height with a detector 7 ft. 
below the floor) are: 

h 


Bui Iding plan area 

Wall thickness 

B q (u,X)* 

r b (u,x)dh 

J o 

0 

20 x 20 

25 

0.45 

0.48 


50 

0.19 

0.20 

o 

X 

o 

25 

0.26 

0.28 


50 

0.086 

0.098 

100 X 100 

25 

0.080 

0.083 


50 

0.018 

0.019 


*More properly, these quantities should be B 0 (u,X)L(d) and / B (i*>,X)l(d)dh, where 

"o ° 

L(d) represents the dose variation with height d above the contaminated plane. If 
this further refinement is included, both quantities are reduced slightly, the integral 
quantity always being reduced more than the multiplication quantity. If the scattering 
wall extends from zero to 10 ft., the integral quantity is reduced by approximately 
8 per cent, while B 0 (w, X)L(d) remains essentially the same. 
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CHAPTER 2 

IN AND DOWN SCATTERED RADIATION - EXPERIMENT 

INTRODUCTION 


Over the past few years several experiments have provided data useful 
in verifying the basement ceiling barrier factor, Bo(Xf,u ), developed in this report^ 
Recent experiments at the Defense Chemical, Biological and Radiation Laboratories 
in Canada, at the Protective Structures Development Center,*“* and at Technical 
Operations Research, ' were performed in such a way as to determine Bg(Xf,u) 
by comparing measured results with and without a basement ceiling slab in position. 
In addition, several other total building experiments have been performed that are 
useful in comparing measured and computed basement data. The latter experiments 
either were not specifically designed for "in and down" measurements or the only 
results available to date are cf a preliminary nature limiting their comparison to 
the total basement protection factor, without isolating the ceiling barrier factor 
component, Bo(Xf,u). Brief comparisons are here presented, however, between 
buiiding attenuation values computed by means of the new barrier curves B^Xpu) 
of this report and the measured results. Of the work cited, only the DCBRL experi¬ 
ment is covered in detail in this report. The remaining experiments have been 
covered in the Conesco report entitled "The Preparation of Simplified Manuals 
for Shielding Analysis."^ 


"In and Down" Scattering Measutcments at the Defense Chemical, Biological and 
Radiation Laboratories 11 5 

The experimental arrangement (see Figure 5) for the "In and Down" 
scattering study at DCBRL^ consisted of measuring the dose distributions produced 
by the radiation transmitted and scattered from a vertical barrier that penetrated an 
adjacent horizontal barrier. A Cesium-137 source collimated to a 10° half angle by 
a lead-filled collimator was utilized as the radiation source. This radiation beam 
impinged on a 3 ft. high by 4 ft. long vertical concrete slab, providing an effective 
ci r cular pattern of incident dose spread over an area 20.5 in. in radius. This area 
was centered at the physical center of the slab on the slab face. Measurements of 
the scattered dose were taken along a vertical traverse positioned 12, 24, and 36 inches 
to the rear of the vertical slab at depths below the top surface of the horizontal slab 
ranging from 3-3/4 to 30 in. Veit'cal slab thickness variations were 0 p>f, 6.45 psf 
plywood, and 34 and 68 psf concrete. The horizontal slabs were of 0, 36.5, and 73 psf 
r> reinforced concrete. The detectors were partially shielded with lead to limit 
their exposure to multiple scattering from the ground. 
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The attenuation attributed to the horizontal (ceiling) barrier was 
determined by taking the ratio of the scattered radiation dose values with the 
horizontal barrier in position to the respective measured values for the zero psf 
horizontal barrier case. Measured attenuation values versus average solid angle 
fraction, w, are shown in Figure 5. This average solid angle fraction, u,for the 
vertical wall subjected to the 20.5 in. radius incident beam was approximately 
determined from u| and <Jg values based on the lower and upper limits of an equivalent 
square area (36 by 36 in.) of incident radiation. Figure 5 also presents the theoretical 
attenuation curves for 36.5 and 75 psf ceilings, as calculated by the method used in 
this report for cesium (.66 mev) radiation impinging perpendicularly on the vertical 
slab. The majority of the experimental points fall below the theoretical curves 
except at small values (w- .175) of <*J, where experimental dose rates were higher 
than those predicted by theory. This increase above the theoreticol curves at small 
u is not expected and cannot be explained at this time. The values of dose upon 
which these data at small u are based are extremely low and any unaccounted air- 
scattered or ground-penetration components even in small quantities could contribute 
significantly at these low rate regions, providing a false reading of this type. Results 
of other experiments discussed later in this chapter did not show this crossover at 
small values of <*. 


"In and Down" Measurements at the Protection Structures Development Center 

Cobalt-60 measurements were performed on a 24 by 36 ft. rectangular 
concrete structure with a basement having an exterior-walI thickness of 50 psf. 

The experiment was conducted both with and without a 50 psf concrete basement 
ceiling slab, so that the basement ceiling effect could be measured directly. In 
Figure 6 is presented the comparisor. between the measured ceiling effect, Bo(Xf,u) 
and that taken from the theoretical values of this report as given in Figure 3. The 
50 psf theoretical curve falls slightly above the 50 psf measured values, giving a 
tatio of experiment to theory of approximately 0.76. The computed values in this case 
therefore were both conservative and close to the 0.815 mean experiment to theory 
ratio determined for c'l the experiments, covered in this report. (See concluding 
section of this chapter.) 


Attenuation of Cobalt-60 Radiation by a Small Iron Cylindrical Structure ' 

A series of measurements were conducted on a 2 ft. diameter cylindrical 
iron structure. The structure was composed essentially of an iron pipe 2 ft. in 
diameter mounted vertically in the ground so that only the upper 2 ft. of its length 
projected above the ground level. The above-ground wall thickness of this structure 
was varied from 1/4 to 1-1/2 in. (approximately 5-60 psf) with floor thicknesses of 
0, 10, 20, and 40 psf iron. Below-ground, i.e., basement measurements were made 
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along the vertical centerline at depths of 1/4, 1/2, 3/4, 1, 2, and 3 ft. Results 
expressed as the ratio of the measured dose rate in the basement with the ceiling in 
place to that measured without the ceiling (zero psf), a direct measure of attenuation 
factor Bo(Xf,u), are presented in Figure 7. 

The measured floor reduction factors vary with solid ongle fraction in a 
manner similar to that of the theoretical curves based on the theory developed in 
this report. The measured values are somewhat lower (more conservative) than 
predicted by theory, as would be expected, since all assumptions and approximations 
used in developing the theory were aimed at producing a conservative attenuation 
value. 

Basement Attenuation — Other Experiments 

A comparison of experimentally measured dose rates with those computed 
by the Engineering Manual^ method (modified for Cobalt-60 concrete), using both 
the floor attenuation data given in that report and those provided by the method 
developed in this report, is presented in Table 4 for all currently reported experiments. 

In general, agreement between expei(mentally measured values and the treatment 
recommended by this study is much better than that provided by the method currently used. 
In particular, the basement ceiling barrier experiments currently being made at NDlJO 
(on a 10 by 20 ft. roofless concrete structure with 46 psf exterior walls and floor slab), 
though still preliminary, provide much better agreement with the method used in 
this report than Engineering Manual values, particularly at the lower depths. The 
same trend holds true for the data in an older study^ ^ of a simple structure with a 
basement, even though those data are somewhat difficult to interpret because of 
inhomogeneity in floor mass. The mean ratio of the experimental results of these 
two full-scale studies to those using the computed method of this report is approximately 
0.70, compared to a mean of 0.815 for all experiments reported. Values calculated 
by the method used in this report are conservative compared with those provided by 
experiment and follow the dose versus depth contours better than previous analyses. 

Protection factor versus depth contour agreement is also good for the 
Kansas State University^ ^ blockhouse end the Technical Operations multistory 1/12 
scale steel modelsJ ' ^ The KSU structure was a square 20 by 20 ft. concrete 
house of 60 psf walls and 56 psf floor and roof; the T/O model was of a 36 by 48 ft. 
six-story windowless building. The absolute measured level of the dose in the 
basement of the KSU blockhouse is higher (prediction less protection) than that 
computed by the method of this report -the only time that this has occurred in a 
full-scale structure. Dose rates in the basement of the 80 psf wall and floor steel model 
(T/O) also were higher than predicted. Reported misalignment of the structure with the 
basement and edge scattering of radiation from the relatively thick basement ceiling 
s'^b that projected above the ground level indicates, however, that the absolute value 
of ! e d ~>se measured is of questionable validity in this case. 
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! TABLE 4 

( Dose Rate in Belov/ Ground Regions from Ground Based 

Sources of Contamination 



>epth, ft 

Experiment 

Engineering 

Manual 

Method of 

This Report 

NDL 10 

1 

.0034 

.0035 

.0052 


2 

.0033 

.0023 



3 


.0016 



4 

.0029 

.0010 

.0043 


5 

.0026 

.00068 

.0035 


6 

.0024 

.00048 

.0028 

Simple structure 

1-3/4 

.0019 

.0010 

.0025 

with basement^ ' 

4-1/4 

.0018 

.00098 

.0023 


6-3/4 

.0017 

.00101 

.0026 


9-1/4 

.0015 

.00040 

.0020 

12 

KSU blockhouse 

1/4 

.0041 

.0026 

.0025 


1-3/4 

.0048 

.0018 

.0025 

Multistory steel 
m „del i3 ' H 

4-3/4 

.0051 

.00082 

.0026 

20 psf walls 

3 

.0195 

.0142 

.0259 

20 psf floors 

10 

.0185 

.0055 

.0185 

15 

.0140 

.0035 

.0106 

20 psf walls 

1 

.00094 

.00164 

.00101 

80 psf floors 

3 

.00115 

.00133 


7 

.0016 

.00075 



10 

.0018 

.00052 

.00176 


15 

.0018 

.00034 

.000127 

80 psf walls 

1 

.00085 

.00072 

.00041 

80 psf floors 

3 

.0030 

.00058 

.00051 

7 

.0033 

.00034 

.00055 


10 

.0034 

.00023 

.00070 


15 

.0026 

.00015 

.00056 

40 psf wails 

3 

.0022 

.0034 

.0040 

50 psf floor: 

6 

.0023 

i .0022 

.0040 

9 

.0021 

.0017 

.0048 


15 

.0022 

.00087 

.0033 


t 
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Summary of Comparisons with Experiment 


The method o: this report for computing basement protection factors, using 
the value Bo(Xf,w) in determining the effect of the basement ceiling barrier, provides 
results that agree much better with measuit’d basement attenuation than does the 
Engineering Manual method, particularly as detector depth below the ceiling slab 
increases. When comparing computed values based on the method used here with 
measured values, nearly all the computed values indicated slightly less attenuation 
than the comparable measured values. The computed or theoretical values therefore 
are conservative as compared with measured values. 

Figure 8 illustrates the comparison between measured and computed values 
of Bo(Xf,u) of this reporr for the three experiments providing data for basement 
positions both with and without the ceiling slab in place (i.e., a direct measurement 
of barrier factor B' 0 ). All values fall on the conservative side of the figure except 
for a few of the DCBRL points token at small solid angle fractions (u <.175), where 
the effects of any radiation component not accounted for (skyshine is suspected) in 
these low dose rate areas would provide a significant part of the measured dose and 
would therefore seem to reflect a major decrease in the floor slab's effectiveness. 
Figure 9 shows the distribution of the ratios of all measured basement attenuation 
experimental values compared with that computed using the theory of this report. 

At ratios greater than 1.10, the tail on this graph can be attributed almostentirely 
to measurements at small solid angle fiactions for the DCBRL experiment mentioned 
above, to the KSU data, and to the questionable 80 psf wall and floor T/O model 
multistory structure data. A statistical analysis of this distribution is shown in 
Figure 10, which clearly indicates the mean value of the experiment to computed 
dose ratio is 0.815 with a sfcndard deviation, a, of 0.135 for all the experimental 
data. 


Since Figure 10 represents a sizable number of experimental data 
(approximately 350 experimental points), it is recommended that the mean ratio of 
experiment to theory of 0.815 be used to reduct the computed barrier function 
(Figures, 2, 3, and 4) of this report in order to provide better agreement. Clearly, 
it is not enough to multiply tire computed values by this ratio, because the attenuation 
of a zero thickness floor must remain 1.0. In Figure 17 of Chapter 3, the floor 
barrier attenuation for the 1.12 hr. fallout spectrum has been corrected by multiplying 
by 0.815 for floor thickness greater than 20 psf, by .91 for floor thickness of 10 psf 
and faired into a vo'ue or 1.0 for floor thicknesses of zero psf. This chart thereby is 
brought into better agreement with expected experimental values. It is further 
recommended that a similar procedure be used for cobalt and cesium. 

A report entitled "Dose Albedo and Transmission Coefficients for Cobalt-60 
and Cesium-137 Gamma Rays Incident on Concrete Slabs' 1 ^ by M. J. Burger and 
E. E. Morris was published when this report was nearing completion. The attenuatior, 
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Figure 8 
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of parallel slant incident radiation for the values used in calculating the floor 
attenuation for this work at 1.25 mev is compared with *hot of Spencer^ ct 1.25 rr.ev 
and Raso at 1.00 and 2.00 mev in Figure 11 for o typical case. This figure 
presents the comparable attenuation curves as a function of the cosine of the slant 
incidence angle 0 for a slab of one mean free path thickness. This latest work of 
Berger and AAorris yields lower values of attenuation (i.e., the transmitted 
radiation level is lower). The Berger end Moms attenuation coefficients ,in general, 
are lower by 10 to 20 per cent than an interpolated 1.25 mev cur/e for Raso’s work, 
the greatest difference occurring in the range 0.3 < cos 0 <0.6. The attenuation 
factor Bg(Xf, u) proposed in this report for use in basement calculations was based 
on the transmission data of Raso. It was further proposed that the Bo(Xf,«) results 
be reduced by o factor of 0-815 to bring theory and experiment into closer agreement. 
If Berger and Morris transmission data hod been used in place of Raso's, the 
theoretical value of Bo(Xf,u) would have agreed nuch more closely with experiment 
without the 0.615 factor. This confirms the necessity of the 0.815 reduction factor 
to bring the analysis in line with experiment. 
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Supplement to TR-20-Vol. 2 


CHAPTER 3 

DRAFT SUPPLEMENT TO TR-20 Vol. 2 and PM100-1 
FOR IN AND DOWN CORRECTION 

IMPROVED BASEMENT CALCULATION 


INTRODUCTION 


Recent studies have indicated that more accurate protection factors may 
be achieved in basement areas if the attenuation of the basement ceiling is con¬ 
sidered to be a function of both ceiling thickness and structure size. In this memorandum 
we shall summarize briefly the method of calculating the attenuation provided by a 
basement ceiling and will provide the charts and graphs necessary to make this 
calculation. 


Floor Barrier Corrections, &XJ(Xf ) 

In the current edition of the equivalent building method, the floor barrier 
effect is represented only as a function of the floor thickness. Recent experiments 
indicate that, in some cases, this interpretation may be somewhat in error. New barrier 
functions providing much better agreement have been computed by assuming single 
Compton scatter energy loss in the structure walls and attenuation similar to that 
experienced by parallel slant incident photons through the floor slab. Floor barrier 
effect thus becomes a function of both floor thickness and structure area. 

Figure ll-l of TR-20 Vol. 2 indicates that the basement Pf charts are based 
on a floor barrier of 1.0 or, equivalently, a mass thickness of zero for the floor above 
the detector This was done so that any corrections would be additive. The effect 
of basement ceiling thickness in attenuating radiation fiom ground-based sources of 
contamination is presented in auxiliary charts, one for each basement area. 

New figures replacing Charts 5 through 9 of TR-20, Vol. 2 accompany this 
supplement. When these charts are used to compute Pf for a bsement location, the 
auxilary charts on the upper poition of the figure must be used to correct for the 
barrier effect of the floor above the detector. The auxiliary charts were derived by 
converting for each building area the attenuation values computed as described above 
to equivalent weight of wall barrier factor B e from Chart 2 of the Engineering Manual 
(TR-20, Vol. 1). 
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Chart 9 


Figure 16 


PROTECTION FACTOR BASEMENT AREA: » 100,000 sq. ft 


32 


















CONES CO consultants in engineering science 


Illustrative Example 

What is the protection factor in the 
basement for a rectangular building 
90 ft. wide by 110 ft. length. 

Floor thickness Xf = 50 psf 

Wall thickness X w = 70 psf 

Roof thi ckness X r = 100 psf 

Solution: 


X 



Building area = 90 x 110' = 9900 2 ~ 10,000 sq. ft.; therefore, Chart 8 

/101 

Equivalent basement area Ag = 9900 = 3,430 sq. ft. 

1 o 2 

Equivalent roof area A' = 9900 = 2250 sq. ft. 


Therefore, from Chart 16 (TR-20 - Vol. 2), equivalent roof thickness ~ 152 psf 
Wall equivalent ceiling thickness AX U (Xf) = 195; Chart 8, top 
Enter Chart 8 with x w = 70 + 195 = 265 psf and X Q = 152 

Read Protection Factor, Pf = 180 
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Supplement to PM 100-1 
TR 20 Vol. 1 


DRAFT 

IMPROVED "IN AND DOWN" CORRECTION 


INTRODUCTION 


Recent studies have indicated that protection factors of higher accuracy 
may be achieved in basement areas if the attenuation of the basement ceiling is 
considered to be a function of both ceiling thickness and structure size. In this 
memorandum v/e shall summarize briefly the method of calculating the attenuation 
provided by a ceiling and will provide a new chart of floor barrier factor for ground 
contributions through the floor above the detector to make this calculation possible. 


Floor BarrierFactor, Bq 

In the current edition of the professional manual entitled "Shelter Design 
and Analysis, TR-20 - Volume 1, Fallout Protection and PM-100-1, Design and 
Review of Structures for Protection from Fallout Gamma Radiation," the attenuating 
effect of the floor slab is represented as a function of only the floor thickness. 

Recent experiments indicate that in some cases this interpretation may be somewhat 
in error. New barrier functions providing much better agreement have been computed 
by assuming single Compton scatter energy loss in the structure walls and attenuation 
similar to that experienced by parallel slant incidence photons through the floor 
slab. Floor attenuation Bq thus becomes a function of both floor thickness and 
structure geometry. 

This new attenuating function, Bo(Xq,w), replacing Bq(X' 0 ), is used in 
the same way as that previously used. The expected dose rate in a structure from 
wal I-scattered sources of radiation is computed for a floor slab of zero thickness 
and then multiplied by the attenuation introduced by the floor. Since the attenuation 
afforded by the floor is a function of the solid angle fraction of the scattering mass 
relative to the protected position of interest, as well as the floor thickness (see 
the figure accompanying the illustrative example), an exact calculation would 
require averaging over the total scattering wal! from Uy to An approximate 
method that provides excellent agreement with this more exact procedure is to 
determine and use a value of attenuation Bo(Xo/ u ) for the arithmetic average of 
the solid angle fraction describing the scattering surface. 
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As in the illustrative example, the attenuation introduced by the floor would 
be equivalent to B^(X' q , u ), where: 




1. Functional equations 

C q = negligible because X Q is greater than 350 psf 
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2. Determining solid angle fractions and parameters 








Chart 

3* 

Chart 

5* 

Chart 

5* 


W 

r~ 

z 

e 


u 

G a M 

G(u) 

u 

u 

90 

no 

24 

.82 

. 

.44 



.33 

w' 

u 

90 

110 

10 

.82 

.18 

.82 

m 



u + u* 

— _ u u .60 + .82 

Y~ ^- 


B' (X ,u 
o c 

= 8^(50,0.71) 

= 0.029 

Accompanying chart (Fig. 17) 

B (X h) 
e e 

= B e (70,3) 

= 0.18 

Chart No. 2* 

S (X ) 

w e 

= S (70) 

w 

= 0.66 

Chart No. 7 * 

E(e) 

- £(.82) 

= 1.41 

Chart No. 8* 


♦Chart of PM 100-1 


3. Solution: 


C = (0.18) (0.02?) { T(.079 - .051) (1-.66) + 


.33- .19 


.661 


1.41 


C = .00073 

g 

K = c + C = .00073 +.00000 = .00073 
r g o 

Pp = l/Rp = 1/.00073 = 1370 say 1300 answer 
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Illustrative Example 2 - Ref. Example 5.3, 
TR-20 (Vol. 1) 

Exposed Basement Walls 



L 

100 ft. 

W = 50 ft. 

X 

80 psf 

(first story) 

e 

X ' = 

100 psf 

(basement) 

e 


x » = 

60 psf 


0 

X 

100 psf 



r 
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Chart 

3 

Chart 

5 

Chart 

5 

u 

w 

L 

z 

e 

n 

U) 

G 

s 

G 

a 

u” 

u 

50 

100 

23 

.50 

.46 

.47 

.39 

.036 

u 1 

u 

50 

100 

11 

.50 

.22 

.70 

.28 

.069 

u 

u 

50 

100 

3 

.50 

.06 

.91 

.11 

.029 


Soluti on: 


a. 


u 


u‘ +u" 

u u 
- 


.70 + .47 

2 


.59 


B (X ,H) = .095 
e e 

S (X ) = .73 

w e 

E(e) =1.34 


C 

9 


G (u 1 ) 
a u 


G (u ) 
a g 


S (X ) 

w e 


+ 


G (o' ) - G (u ) 
s u s u 


i (X )E(e)> B (X ,H) 

we i e e 


C 

9 


(.069 - .029) (1-.73) + (.28 -.11) (.73)0.34) 


C =. .0168 
9 


.095 
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c. Roof Contribution 

v = X + X’ = 160 psf 

o r o 

i 

C (u\X ) = .0039 

o U 0 

P f = 47 

Detector heights above ground are treated in a similar manner 
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CHAPTER 4 

FINITE FIELDS OF CONTAMINATION 


INTRODUCTION 


Many shelter areas in this country depend in part for their worth on the 
fact that they are subjected to limited rather than infinite fields of contamination. 

A common situation, that of row housing in urban areas, for example, would 
subject the sheltering structures to radiation originating from a finite field of 
contamination perhaps 50 ft. wide on both sides of the structure. These structures 
would then be as much as three times more effective than an identical structure 
surrounded by a full field of contamination. Parker^ estimates that fully half 
of the basement shelter areas surveyed at the time his report was published were 
subjected to field widths of less than 60 ft. 

Two areas of inquiry arise when we evaluate the methods currently used 
in predicting the effects of finite fields of contamination. We must first determine 
if the more detailed method of calculation as represented by the "engineering 
method" adequately portrays the physical situation. Second, we must ask whether 
or not the simplified "equivalent building method" agrees with results obtained 
by means of the more complex method. The agreement between the two methods in 
their present form has been examined^ previously and recommendations f or improving 
upon this agreement have been made. The required changes have been held in 
abeyance until an adequate answer to the first area of inquiry (the accuracy of 
the more complex method of calculation) has been determined. Our purpose in 
this chapter is to examine this accuracy by comparing if with the best information 
available. 

There are few experimental data with which to evaluate the methods 
proposed for computing the shelter afforded against finite fields. Several years 
ago, an experiment was conducted to determine the effect of finite fields of 
contamination on above-ground locations. This experiment used the modeling 
technique, with iron as the attenuating medium, and simulated contamination in the 
form of rectangular areas. Eisenhauer* analyzed these data in detail, concluding 
that "comparison of experimental and calculated results shows good agreement for 
both structures except for the first story." A second, brief experimental series, con¬ 
ducted on full-scale concrete structures by McDonnell, using circular annuli, 
again indicated fair agreement between experiment and theory. Both experiments 
were exclusively for upper-story locations, where both direct and scattered radiation 
affected the dose measured in approximately equal amounts. Since neither of these 
experimental series was rigorous in the examination of finite fields of contamination 
and because so many shelter areas seem to be affected only by finite fields of 
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contamination, a further evaluation was undertaken. Eisenhauer recently com¬ 
pleted a computer compilation program, using Berger and Morris' Monte Carlo 
transmission coefficients' ^ for 1.25 mev prediction. This program is capable of 
computing the raaiction originating from any arbitrarily located square patch of 
contamination through any rectangular wall area to an arbitrarily located detector. 
The program, therefore may be used to generate numerical "experiments" with 
which to compare the existing calculational method. 


The Engineering Method of Shelter Computations from Finite Contamination Fields 


The engineering method of shelter computation described in the manual 
"Design and Review of Structures for Protection from Fallout Gamma Radiation"^ 
presents rules for computing the ground-based dose contributions for structures 
subjected to radiation from both infinite and limited fields of contamination. In this 
method the radiation arriving at a point within a structure is subdivided into three 
components: radiation that (1) has passed directly through the building walls without 
scattering, (2) has been scattered by the wall', and (3) has been scattered by the 
atmosphere. Non-wall scattered or direct radiation from an infinite field of con¬ 
tamination is determined by multiplying the cumulative angular distribution of 
non-wall scattered radiation, Gj^h) 1 as viewed from the point of interest in 
the structure by the height-dependent wall barrier factor B(X e ,h)i and by the 
fraction of ladiation not scattered by the structure walls, (1-S W ). Calculations 
for the finite field case are similar except that the cumulative angular distribution 
of non-wall scattered radiation must be differenced to account for the finite field. 

The wal I -scattered component for an infinite field case is determined by 
multiplying the cumulative angular distribution of radiation scattered from the 
structure walls, G s (u) , by the height-dependent wall barrier factors B(X e ,h) 
the shape correction facror, (E), the fraction of emergent radiation scattered, (S w ), 
and a floor barrier factor, B 0 (Xf) or ceiling factor, Bo(X c ) , where applicable. 

The floor barrier factor is used for attenuating wal I-scattered radiation reaching 
the detector from the floor below in a multistory structure whereas similarly, the 
ceiling barrier factor is used for wall-scattered radiation reaching the detector 
from the story above For a finite field case this procedure is modified by applying 
a different wall barrier factor, Bu^u^Xg) , which is expressed as a function of 
the solid angle of the field of contamiation as viewed from the structure's wall at 
midsto'y height in place of the infinite-field wall-barrier factor, B(X e ,h). 

The atmospheric scattered component for both the finite and the infinite 
field of contamination cases are determined by multiplying the cumulative angular 
distribution of skyshine radiation, G a (w), , by the wall barrier, B(X e ,h) , the 
fraction of radiation not scattered in passing through the structure wails, (1-S W ), 
and the ceiling barrier, Bo(Xe) , when appropriate. The assumption that the 
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atmospheric scattered component for a finite field is nearly identical to that of 
the infinite field case should, in practice, be fairly accurate. Though adjacent 
buildings may block the direct view of part of the field of contamination, however, 
the air-scattered radiation physically blocked by the adjacent buildings will be 
approximately balanced by radiation scattering from the walls of these buildings. 
Also, it is expected that in most realistic cases .here will be fallout on the tops 
of the adjacent buildings as well as in the region behind them; hence the 
angular distribution of the radiation that is not blocked by adjacent walls will 
remain the same. 


The three radiation components: direct, wail-scattered, and air-scattered, 
are then summed to give the total dose contribution. The general equation for the 
total ground contamination within a multistory structure in a limited field of 
contamination is presented, together with a sketch of the idealized building 
arrangement, in Figure 18. 


Charts in the Engineering Manual for determining the geometric and 
barrier reduction factors used in structure shielding calculations were derived from 
the basic data on radiation penetration developed by Spencer'^ f or 1. ]2 hr. fallout. 
Chart 7 in the manual, for the fraction of emergent radiation from a wall barrier 
that is scattered, (S w ), is, however, for Cobalt-60 radiation. Since most experi¬ 
mental verifications of the Engineering Manual calculation procedures have been 
and are being performed using Cobalt-60 radiation, and since the NBS computer 
compilation program^ against which we wish to compare the engineering method 
uses the 1.25 mev (Cobalt-60) Monte Carlo transmission coefficients of Berger 
and Morris, the Engineering Manual method has had to be extended to 1.25 mev 
radiation in order to permit accurate comparisons. Cobalt-60 curves for the 
barrier shielding reduction factors, Bo(Xf) Bo(X e ,3‘), and Bq(X c ) are presented 
in Figure 19a. The data in this figure has been calculated in the same way and 
are equivalent to Chart 1 of the Engineering Manual, in which Case 1 applies to 
a horizontal floor barrier, Bo(Xf), Case 2, the exterior vertical wail barrier, 

B(X e ,3'), and Case 3, an overhead or ceiling horizontal barrier, Bq(X c ). Figure 19b 
presents the Cobalt-60 overhead or ceiling horizontal barrier Bo(X c ). Figure 19b 
presents the Cobalt-60 exterior wall attenuation effects for detector heights to 
300 ft., and is equivalent to the Chart 2 fallout curves of the Engineering Manual. 
Cobalt-60 directional responses from ground sources of contamination Gj, G s , 
and G a are given in Figure 20a, which has been calculated in the same way and 
therefore is equivalent to Chart 5 of the Manual for 1.25 mev radiation. Figure 20b 
similarly extends the directional response for direct radiation G(J to detector heights 
of up to 100 ft., and is equivalent to Chart 6 of the Manual. 


Barrier shielding effects B^Xf), B e (X e , and Bq(X c ) for Cobalt-60, os 
given in Figures 19a and 19b, were developed directly from Spencer's^* L(X), 
2W(X,d), and S'(X) curves respectively, and agree closely with those developed by 
LeDoux.^l 
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Figures 20o and 20b for directional response Gj(u,h), G s (u) ond G Q (X) 
were developed from the basic work of Spencer. *6 LeDoux^l also has developed 
Cobalt-60 curves for the directional responses, which in general agree well with 
the curves of Figure 20 except for values of G<j at detector heights above 3 ft. 

Values for the direct angular distribution G^(u,h) in Figures 20a and 
20b were developed from the equation: 


cos 9v 

r 

j (d,cos9)d(cos0) 

G>,h)= 9 -,- 

I K (d,cos9)d(cos9) 


1 o .1 

(d,cos9)d(cos9) - J<(d,co$0)d(cos9) -J^(d,cos0)d(cos0) 
-1 -1 co$9x, 

^,, , (d / cos9)d(cos9) 


which in Spencer's^ notation is identically equal to 


G d (u,h) = 1 


5(d) 

Lid) 


- L (d,u) 
a 


where t' (d, cos 9) is the angular distribution and S(d), L(d), and L a (d,u) for Cobalt-60 
are taken from Spencer's^ work. S(d) divided by L(d) is the ratio of skyshine to 
total (adiation at a distance d above a uniform infinite source field. G s (cobalt) 
was calculated from the equation used for G s (fallout), using Sa(d,u) curve for 
Cobalt-60 at d = 3 ft. Thus: 


G = 0.5 (1-S (d,u) 
s a 


Similarly, the cumulative skyshine contribution 


G = 0.168 G (u) 1 + (1/2 - G (u) 


where 0.168 is twice the skyshine fraction of the total radiation at 3 ft. altitude. 
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Agreement for Gjfu^ 1 ) between values calculated in this manner and 
curves developed by LeDoux are excellent up to u ; .90. For values of o -.90, 
the values of Figure 20 are somewhat higher than those of LeDoux (for example, 
at u - .99, Gd (Figure 20) .14, G^(LeDoux) - .105). Similarly, agreement for 

the cumulative angular distribution of wall-scattered radiation G $ between 
Figure 20 and values previously derived by LeDoux is excellent, Conesco values 
being slightly higher at o .75 (for example, at u = .99, G s (Figure 20) = .013, G s 
(LeDoux) = .01 1). The skyshine cumulo ve angular distribution G a also shows 
good agreement between the values shown in Figure 20 and those previously 
derived by LeDoux, with values of Figure 20 slightly less for u .75 and slightly 
higher for u .75. Figure 20b presents the direct response Gj for heights to 
100 ft. and agrees fairly well with values developed by LeCoux. Agreement at 
Gd(u, 3') as discussed above, was good. Agreement at heights greater than 3 ft, 
however, is not as good, particularly at small values of o. This is illustrated in 
Figure 20 by including both curves at Gd(u,h) equal to 0.50 and 0.75. 

Engineering Manual Chart 9, used for obtaining the finite field wall 
barrier factor Bu s ( u s ,X e ) is for 1.12 hr, fallout. Fallout values for Bu s ( U s/X e ) 
were used in this study for calculating 1.25 mev structure reduction factors for 
comparison with the computer Monte Carlo results. This was done for two reasons: 

(1) a Cobalt-60 barrier chart for limited strips of contamination had not been 
developed at the time this study was conducted; (2) the error introduced by using 
the fallout chart values in place of Cobalt-60 should be small at the mass thickness 
(9, 36, and 144 psf) used in this study. The error in the finite wall-barrier factor 
in using failoui values instead of Cobalt-60 can be approximated by considering 
the differences between fallout and Cobait-60 infinite field wall barier factors. 

For the infinite field case the Cobalr-60 barrier factor is about 5 per cent less than 
the fallout value at 100 psf, and 25 per cent at 200 psf. Since the wall thickness 
values used in tfvs study are 36 psf or less, except for one point at 144 psf, errors 
introduced are expected to be small. 

Computations of the Engineering Manual type were performed using data 
for 1.25 mev on a series of simple structures for comparison with corresponding 
values determined by the NBS compilation prog _ om for finite fields. T.ie series 
v/as based on a single-story structure of the blockhouse type. Building floor area, 
widtn-to-length ratios (eccentricity), and wall thickness were varied to provide 
data for a comparison with the compilation program over a wide rangr* of structure 
parameters. The width of the rectangular field of contamination surrounding the 
structure was varied from a narrow strip only 10 ft. wide to a maximum width of 
90 ft. Gne group of structures provided buildi ig iength-to-width ratios from 1 to co. 
The bases of these buildings v ere 20 by 20 ft., 20 by 40 ft., and 20 ft. by co, 
calculations being made for each for three wall mass thicknesses (9, 36, and 144 psf). 
A second group consisted of two square buildings 20 by 20 ft., and 40 by 40 ft., again 
computed for wall thickness values of 9, 36, and 144 psf. The building height was 
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10 ft. ond the detector was centrally loccted at 5 ft. height for all cases. Computed 
values of the direct ond scattered radiation for the two groups of structures are 
summcrized in Tables 5 and 6. 

In certain cases the values of direct radiation presented in Tables 5 and 6 
are slightly different from those which wouid be computed using the data of 
Figures 19 and 20 because the values of Gq used in computing these results varied 
slightly from the refined values presented in Figure 20. The differences in direct 
radiation caused were generally less than 10 per cent and do not affect the 
conclusions presented in this chapter. 

Shielding calculation of the Engineering Manual type as illustrated in 
Figure 18, for this series of single-story structures subjected to finite rectangular 
fields of contamination may be simplified to the following equation: 


Dose = direct radiation component + scattered radiation component 


D 

o 


B(X ,h) 
e 


G d (o. ,h) - G d ( U *,h)'j (1-SJ ♦ B os ( Us ,X e ) [g (», + G < V > 


S E 

W 


The skyshine component is omitted from this equation so as to permit direct 
comparison with the NBS Monte Carlo compilation method. The NBS method 
neglects the effect of the atmosphere and assumes that all radiation enters through 
the detector story wails of the structure (i.e., there are no radiation components 
from the story below or above). The direct radiation attenuation factor, B(X e ,5‘), 
was taken from Figure 19b, and the angular distribution, Gj(u ,5') and G<j(u*, 5'), 
generally from Figure 20b. Two complications, however, enter into the calculation 
of the direct radiation components as the building width/iength ratio approaches 
zero. The first arises in calculating solid angle fractions for an infinitely long 
building and *he second is introduced by the tacitly assumed equivalence of circular 
and rectangular solid angles when com, ing angular distributions. For an infinitely 
long structure, the solid angle fraction <j equals _ tan “I ( W ) r where W is the 

* 2Z 

structure width and Z is the detector height. As for the second complication, one 
finds that in differencing G^O-^h) for a series of buildings with fixed width and 
varying length and a fixed finite field width, the value of A Gd first decreases 
with increasing building length as one wouid expect, then as the building increases 
further in length, A Gd increases. As building length increases toward infinity 
the direct contribution and hence A Gj should asymptotically approcch a minimum 
at a length-tc-width ratio of infinity and should not increase. In the actual 
calculations for direct radiation for a given source field and building width, AGd 
was plotted versus building eccentricity until the minimum was reached. Values 
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beyond this minimum were extrapolated with the extrapolated curve becoming 
tangent to o constant value as V-/L - 0. Thc^. extrapolated values were as much 
as 10 per cent less than hie calculated — value. 

The NSb Compilation Prog^cm for finite Fields of Contamination 

2q 

Eisenhas/er has recently completed a program to compute the effect of 
a finite field of contamination using transmission and scattering coefficients com¬ 
puted by the Monte Carla method by Berger.^ The advantages of this type of 
calculation are, first, that it provide- a completely different method of calculation 
starting with different assumptions and a different set of parameters and input data, 
against which the current limited stri>; calculation may be compared and, second, 
since this program utilizes transmission coefficients computed by the Monte Carlo 
method it may truly be thought of as a number! cal representation of the cctual 
experiment. This code has been designed to c mpute the dose contribution to a 
detector located behind a vertical wall, from contamination locc‘ * in a rectangular 
area, calk 1 o source pofh (with its axis parallel and perpendicular to the wall) 
through i verficu' rectangular sec den .tie wall. Detector, wail section, and 

source patch may be iocatea at any arbitrary position in relation to each other, 
individual cont. but?on-, from each source patch representing differential areas of 
a field are computed a;,d presented separately. Since the calculation is designed 
for finite fields of contamination, air scattering and attenuation are neglected. 

That this effect is negligible for finite fields of the size of interest may be quickly 
demonstrated. The dcso late from a limited circular field of co^tamination to a 
detector located a height h above it is: 


ft 

o 

D * 2ttS q dp 

° P 

h 


where p 



toral cross section ior cir- J/450 ft. S.T.P. 

specific irradiance of the source (raas/hr curie) at 1 ft. 

source density (curies/sq. ft.) 

detector height (it.) 

outer siant radius of the field (ft.) 
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j 


vAieh in t is absence of scattering {all collisions are assumed absorption, clearly 
placing an upper limit on the effect of the atmosphere) reduces f? 


D = 2a 


S o q 


— 

E 


E,G»H)- 


1 


(pR )! 

°J 


where 


E 


1 


the exponential integral of the first kind 


if the ■c;««osoh*re is totally neglected p is zero and the corresponding equation becomes: 


D = 2. S o qln ^ 


The maximum error introduced by neglecting the atmosphere would thus become 


Error = 


E,(ph) - £,(pR o ) 
- 

In £ 


Field Radius (ft.) = 0 
Error (%) = 0 


10 20 30 

1.2 2.2 3.0 


50 70 100 

4.3 5.5 6.8 


The dose resulting from uncollided photons is computed in Eisenhauer's 
program by first determining if a straight line, from the center of the source patch ro 
the detector, passes through the wall area (see Figure 21). If it does, it is assumed 
that all the contamination in the rectangular source patch is concentrated at the 
center of that ^tch. The dose at the detector is computed by reducing the output of 
that contamination by the distance to the detector squared and exponential attenuation 
along the slant distance through the wall (determined by the straight line connecting 
the source patch to the detector location). Scattered radiation is computed in a 
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of the rec’mqirici wall crea of interest and then drewing a second line from this 
point !o the defector 'location. The o.'gle mode by these two lines is then 
determined one a transmission coefficient for scattered radiation is interpolated 
from the TAonfe CsHo transmission foci's of &erc?r r The radiation emanating from 
ihe scarce patch is then reduced by the radius to me cer-ter c? th= •*«!! ore- squared, 
the transmission vOeftTci c r -t o°d the radius to the detector square". 


Clearly, both the waii oris score-? patch grid *?«»* be ct-:ose,? with cure so 
that furih-j.: in grid size will yield a negligible charge in resists, 

fxpipraicry a. eolations indicate that for a detecsc-i one or more heights behind 
the wall, locofed at the wall •mid-height, a square wo'.* grid -,>ne-hc{* I height 
in dimension together with source patches of a similar size for the firsi l «c wail 
h .rights av/cy from the wcil ana twice this size beyond, aroduces results that are 
within i per cent of those produced b, a grid one-fourth as large. Eisenhauer^ 
estimates that this program produces results which are accurate to about 5 per cent if 
the grid sizes are properly chosen. 


After preliminary evaluation of the program the above-mentioned grid size 
was selected as adequate for these studies. The resultcnt dose contribution from o 
source pctch on the right of the centerline through the de'ector perpendicular to the 
wall are: 

1. The direct contribution — the dose contribution that proceeds 
directly to the de’.ector without scattering. Since the detector 
is located at th*. vail mid-height, contribution arises through 
the lower half of the wall only (direct). 

2. 7 ie <vull-scatter contribution through a wall slab in the lower 
hnlf of the wall to the right of the centerline (wai! scatter lower). 

3. The .ymmetric wall-scatter contribution through a wall slab 
symmetrically located on the left of the centerline (symmetric 
scatter low-.-'). 

4. The wall-scc'rer contribution through a wall slab in the upper 
half of the wall to the right of the centerline (wall scatter upper). 

5. The symmetric wall scatter contribution through the upper wall 
slab symmetrically located from (4) to the left of the centerline 
(symmetric scatter upper). 

Although the balance between direct end scattered radiation is set by 
the wall thickness, the division between the modes of scatter penetration seems 
quite general (at least for the cases investigated, between one-quarter and four 
mean free path thickness). For source patche r located close to the wall the wall- 
scatter from '.he upper half of the wall is ger.eially from 5 to 20 per cent of wall 
scatter from the lowc' portion. As the distance from the source patch to the wall 
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increases, this r'ivisicn becomes more neai'y equcs, «iih values generally within 
10 per . snt of one another at 80 ft- distance. Scattered radiation arriving through 
symmet..cclly located dctecto:s generally behoves similarly when toe source patch 
end wail slob ere on er near the centerline. As either the slab or the source patch 
moves off-center, hc---.sever. ;;>> symmetric scatter radiation contribution decreases 
mere markedly. Results obteint or a 20 by 20 ft. square building with a 5 ft. 
detector height are illustrated in figure 22 for each radiation component. The 
important feature ct Figure 22 is the s'^ope of the curves cs the width of the 
contaminated area increases. Ecch component of transmitted radiation increases 
■*i almost the same way. 

The NBS Monte Carlo compilation program is both extremely flexible 
and easy to use (though summarizing the results in a form usable for this study 
necessitated much manual 'abor in generating the required number of "numerical 
experiments” w:*h which the cresent method c-f calculation could be evaluated. 

As an initial stud-, only the single detector height or 5 ft. was chosen, together 
with c- 'taicfure woi'i height of 10 ft. Building size and wall thickness were then 
varied and the data zed :w t’-mt the finite rectangular limited fir kr 

of contamination surrounding a structure could be cc'culo^d. The data obtainc-d 
from this series of calculations are summarized in Table 7. 


The two methods of calculation compared 

As illustrated earlier in tr>*s chapter, rhe maximum effect of air scatter 
and attenuation on the radiation originating from contaminated areas up to iOO ft. 
wide is only approximately 7 per cent. If both methods of calculation produce 
the same results, agreement should be within this 7 per cent, tbs- engineering 
manual results being lower rr.cn those computed by rhe NBS Monte Carlo compilation 
method. Figures 23, 24, and 25 present this comparison for total dose a function 
of width of the contaminated field. Figure 22 represent', a cose in which the 
eccentricity (W/L) of the building is varied from zero to one whh all other 
parameters held constant. Figure 24 presents a series of caint.-lotions in w'-'icl-. on.y 
the wall thickness is varied; Figure 25 represents two cases ider.ficol except for 
building plan area. These graphs indicate excellent agreement in iota! dose os 
paedicted by the two methods. The mean value of the ratio of the total dose as 
computed by the NBS methods to that computed by the Engineering Manual method 
I for all cases investigated is .98 (when *he ratio is computed at V'^h = 2,3,4,6,8, 

10,12, 14 for each case investigated), with average values rancyng from .91 to 
1.04 for any one structure design, i? was expected that this ic'io would be slightly 
greater than one (because atmospheric effect: were neglected in the NBS method, 
rather than less than one. This minor difference mcbably wa f caused by the accuracy 
to which the engineering method graphs may be read. Since the ratio is so cl se 
to one and since the structure investigated was varied in size, wall thickness, and 
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Normalised so as to make the response 3 ft. above or infinite field of contominati 
1.0 if *he effects of the atmosphere were included. 
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eccentricity, it may be safely concluded therefore that at first-floor locations as 
typified by detector height of 5 ft., the current method of computing the effect 
of finite rectangular fields of contamination is adequate. 

Unfortunately there was not enough time to generate a similar group of 
Jata for comparison at other detector heights (for example, in bcsement areas and 
upper-floor locations) in order to make the conclusions more general. In an attempt 
to estimate the effects at other detector locations, the scatter and direct dose 
compared component by component. A surprising result was uncovered: although 
the ratio of either direct or scattered dose as computed by the NB5 Monte Carlo 
Method to that computed by the engineering method was very nearly constant as 
a function of field size (significant variations from the mean value occurred only 
within the first 15 ft. of field size), this ratio varied significantly from one. 

Table 8 presents a summary of these data computed, as before, for each case investi¬ 
gated at field widths of 5, 10, 15, 20, 30, 40, 50, and 70 ft. 


TABLE 8 i 

Mean ratios of scatter and direct radiation as computed 
by the engineering method (10 ft. structure height, 

5 ft. detector height) 


5 ! 


_ 

rati 

os 


W 

Parameter 

varied 

Structure 

(ft.) 

Wall thickness 

(psf) 

Direct 

Scat. 

Old 

Indicated 


20 x 20 

9 

1.13 

.513 

.178 

.11 

Wall thickness 

20 x 20 

36 

1.45 

.751 

.475 

.34 


20 x 20 

144 

1.92 

.846 

.828 

.70 


20 x 20 

36 

1.45 

.751 

.475 

.34 

Eccentricity 

20 x 40 

36 

1.39 

.714 

.475 

.33 


20 x oo 

36 

1.08 

.800 

.475 

.37 


20 x 20 

36 

1.45 

.751 

.475 

.34 

Size 

40 x 40 

36 

1.48 

.541 

.475 

.27 


Sizable variations between the ratio of direct and scattered radiation occur 
as a function of wall thickness, eccentricity, and structure size. Because data 
available are limited, it is difficult either to pin down the reasons for this variation 
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or to recommend a change. Since the ratios of total radiation computed by both 
methods were so nearly one, however, the ratio of scattered to total radiation, 
S w (X e ), used in the current methodology becomes suspect. The values of this 
parameter required to make the ratios of Table 8 equal to one are presented in that 
table and in Figure 26, together with the values currently used. Because of this 
discrepancy we must therefore caution that at this time the conclusion that the 
engineering method of computation is valid for limited fields of contamination may 
be made only for locations on the first floor of a structure . Other locations, where 
the relative importance of direct or scattered radiation is considerably different 
(as in basement areas or upper-floor locations where the field is shadowed by a 
thick floor) must be investigated further. 
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FRACTION OF TOTAL RADIATION THAT IS SCATTERED 
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CHAPTER 5 

CONCLUSIONS AND RECOMMENDATIONS 

In this study we have examined the new method ^proposed for computing 
the attenuation of "in and down" scattered radiation and the current method of 
computing the effects of finite fields of contamination so that we might recommend 
improvements and up-dating for OCD documents OCD-TR-20, Volumes 1 and 2, 
"Shelter Design and Analysis." 


in and Down Scattered Radiation 


The new barrier factor B' 0 (X / u ) computed for in and down scattered radiation, 
based on the assumption of Compton first collision scatter in the structure walls, has 
been found to provide significantly better agreement with existing experimental 
measurements. Comparing this factor with approximately 350 experimental measure¬ 
ments, we found that the new method was generally conservative, predicting dose 
rates from ground-based sources of contamination 18.5 per cent higher than those 
measured, with a standard deviation of 13.5 per cent. 

It is recommended that this new barrier factor B'o (X,u ), reduced by 
18-1/2 per cent, replace the currently utilized Bq(X) to improve the accuracy of 
the calculation method. In Chapter 3 of this report we provide detailed recommendations 
on the steps required to effect this change in both the Engineering Manual, 

OCD TR-20, Vol. 1, and the equivalent building method, OCD-TR-20, Voi. 2. 

Finite Fields of Contamination 

A preliminary investigation to determine the effects of finite fields of 
contamination for detectors on "first-floor" locations indicates excellent agreement 
between total doses when calculated by the two different methods, 'vloior variations 
in the scatter to direct ratio were obtained by the two methods. There is, therefore, 
serious doubt as to the accuracy of the recommended procedures for positions where 
the scatter to direct ratio is significantly different from that investigated (e.g., in 
basement or upper-story locations where the detector is shadowed by a heavy floor 
slab). It is recommended that the effects of finite fields of contamination in areas 
where direct radiation does not provide a major dose component be studied further 
in order to develop a simplified method of calculation that will be accurate in all 
locations. 
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General 


The current program has reviewed most of the experimental and 
analytical efforts devoted to shielaVtg from contamination in a "fallout geometry." 
It has become evident during this review that, in addition to the two areas of 
inquiry that are covered in this report, several additional areas should be investi¬ 
gated. These are: the ratio of scattered to direct radiation, particularly as 
applied to finite fields of contamination in below-ground and upper-story shelter 
locations; a preliminary investigation of simplified methods for computing the 
effects of inhomogeneous floor and wall barriers; and a preliminary study of the 
effects of roof wash down (or other decontamination methods) and how these 
effects may be simply expressed (perhaps as an additional equivalent roof thickness) 
in each of the present methods of caltulation. The latter program is suggested so 
that this procedure (of decontamination) may become useful in analyzing new 
building construction when other slanting or design techniques cannot be used. 
Additional effort is required in each of these areas to ensure that the appropriate 
methods of analysis currently used or proposed are accurate, simple, and easily 
usable so that they may be included with confidence in OCD publications. 
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